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ABSTRACT 

We analyze a technique of obtaining turbulence power spectrum using spectral line data along 
the velocity coordinate, which we refer to as Velocity Coordinate Spectrum (VCS). We formalize 
geometrical aspects of observation through a single factor, "geometric term". We find that all variety 
of particular observational configurations can be described using correspondent variants of this term, 
which we explicitly calculate. This allows us to obtain asymptotics for both parallel lines of sight and 
crossing lines of sight. The latter case is especially important for studies of turbulence within diffuse 
ISM in Milky Way. For verification of our results, we use direct calculation of VCS spectra, while the 
numerical simulations are presented in a companion paper. 
Subject headings: 



1. INTRODUCTION 

Turbulence is an essential property of interstel- 
lar medium, which controls a number of key as- 
trophysical processes, including star formation, mass 
and heat transp ort an d propagation of cosmic rays 
(see iSchlickeise r 1999, Va zguez-Semadeni et al.l 120001 
iNaravan k Medvedev .20011. ICho et al.ll2003l . see review 
bv lCho. Lazarian fc Vishniadl2002D . Observational stud- 
ies of turbulence can put constraints on theoretical ex- 
pectations for such processes and clarify their physical 
descriptions. 

The advantage of obtaining the spectra of turbulence 
using the observations of Doppler broadened lines was 
obvious from the very beginning of the research in the 
field (see lLazaria i ([2006) for a recent review) . The prob- 
lem happened to be extremely difficult, however. 

Fortunately, we have seen some progress in this field 
recently. For instance, the VCA techniques for the 
analysis of turbulence from spectral data, developed by 
iLazarian fc Pogosvanl ([2000) (henceforth LPOO), is based 
on the analysis of fluctuations in channel maps. Those 
can be viewed as slices of a Position-Position- Velocity 
(PPV) cube. However, there exists a complimentary 
technique, employing statistics along velocity coordinate 
which has also been briefly discussed in that paper. This 
technique was introduced in the same (LPOO) paper, but 
was not appreciated enough at the time of its introduc- 
tion. In fact, it happened that studying the spectra 
of channel maps was an ad opted procedure well before 
ILazarian fc Pogosvanl ()2000f ) posed a question what these 
fluctuations may mean in terms of the underlying turbu- 
lence. On the contrary, has been studying the fluctu- 
ations of PPV emissivi ty along the velocity c oordinate 
until very recently (see IChepurnov et al.ll2007f l. As the 
result, the introduction of the technique was delayed. A 
detailed study of the technique, that was termed Veloc- 
ity Coordinate Spectrum (VCS) (Lazarian 2004), includ- 
ing the possible effec t of ab sorption were presented in 
ILazarian fc Pogosvanl ()2006f ) (henceforth LP06). How- 
ever, there only a distant cloud is considered there as a 
possible turbulent volume for studies, which exclude the 
situation of not parallel lines of sight, which is typical, 
e.g. for studies of high latitude HI. 



In the present paper we present a formalism that in a 
concise form allows to reproduce the earlier results in the 
field, as well as to move the field forward. The gist of our 
approach is formalization of the role of observational ge- 
ometry, which allows us to write out the basic equations 
of both the VCA and the VCS in a very similar way. 
Our approach also enables us to account for a limited 
resolution, various beam geometries and, more impor- 
tantly, account for both the crossing or parallel lines of 
sight geometry of observations. In other words, unlike 
LP06 we can both consider turbulence studies within a 
distant cloud with size much larger than the distance to 
the cloud and the in volume that surrounds the observer. 
Moreover, our general formalism can easily be extended 
for the intermediate cases, when lines of sight are not 
crossing, but cannot be considered parallel either. 

We note, that the VCS is special among the techniques 
of turbulence stud ies. Compared to the VCA or Veloc- 
ity Centroids (see lEsauivel fc LazarianI (|2005f) and ref- 
erences therein) it uses exclusively data along velocity 
coordinate. In this case the spatial information is used 
for improving statistics only. Our simulations in the com- 
panion paper confirm, that a very few independent mea- 
surements (e.g. ^ 10) is sufficient to get the underly- 
ing velocity turbulence spectrum with enough accuracy. 
This tool is irreplaceable when we do not resolve the ob- 
ject enough to analyze power spectrum in the pictorial 
plane. In addition, our present study makes the VCS the 
only tool accounting for convergence of lines of sight. 

In what follows in Sect. [2] we introduce the basic terms 
of statistics of turbulence, in Sect. |3| we describe the 
signal at a spectrometer output in terms of underlying 
velocity and emissivity fields, in Sect. |3|we consider the 
basic statistics common for VCA and VCS, the power 
spectrum along the velocity coordinate is presented in 
Sect. [5] and in Sect. |7|we consider its asymptotics. The 
variants of the geometric term are presented in Sect. [6l 
Some needed here facts regarding velocity and density 
statistics are presented in Appendices |A] and |b1 the in- 
fluence of the regular velocity shear is considered in Ap- 
pendix [C] We discuss our findings in Sect. [H 

2. UNDERLYING VELOCITY AND DENSITY 
TURBULENCE 
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We assume here, that turbulent velocity and density 
are homogeneous and isotropic and admit a statistical 
description in terms of correlation (or structure) func- 
tions and corresponding power spectra. The latter pair of 
functions is related to each other through Fourier trans- 
form. 

For a vector field like velocity, correlation function and 
correspondent power spectrum take a tensor form. Power 
spectrum is different for solenoidal and potential field 
components. For each case spectrum's tensor properties 
can be represented by a factor, depending only on angles. 
Another, scalar, factor is responsible for energy distribu- 
tion in the wavenumber space. 

Because of the self-similarity of turbulence, we usu- 
ally expect a power law dependence for power spectrum. 
This power law can be either "long wave dominated" (or 
"steep") if the respective spectral index is greater than 
3, or "short wave dominated" (or "shallow") otherwise 
(see LPOO). To keep total energy limited, both types of 
spectrum must have a cutoff: at low wavenumbers for the 
steep one, and at high wavenumbers for the shallow one. 
We assume the velocity spectrum to be steep, as the shal- 
low velocity spectrum is not physically motivated. On 
the contrary, for density spectrum we consider both steep 
and shallow spectra. The shallow spectrum emerges for 
high Mach numbers (see Esquivcl & Lazarian 2005). 

The quantities we deal with in spectral line observa- 
tions are velocity and emissivity. The latter can be 
proportional to density (like in emission lines of neu- 
trals) or density squared (dielectron recombination lines 
in plasma). The regime of squared density modifies the 
underlying density spectrum as the emissivity spectrum 
is an auto-convolution of the density one in this case: the 
asymptotic slope of the emissivity spectrum is the same 
for steep density spectrum and different, = 2a p — 3, 
for the shallow one (see Appx. IB.1|) . 

We are interested in velocity statistics, and emissivity 
is considered only because it affects the quantities, de- 
rived from observational data, that we use here. This 
infiuence is less for steep emissivity spectrum, and could 
be ignored (apart from the normalization), but for the 
shallow spectrum it should be accounted for. 

Because of non-linear character of the statistics we use 
Eq. [21], we can not directly recover the velocity spec- 
trum. But if it is parametrized with the spectral index, 
spectrum amplitude, and cutoff wavenumber (related to 
the injection scale), we in principle can recover these pa- 
rameters. 

3. SPECTRAL LINE SIGNAL AT THE SPECTROMETER 
OUTPUT 

The presentation of the techniques in both LPOO and 
LP06 starts from presenting the underlying turbulence 
statistics and proceeds with introducing the observables. 
In such presentation of the VGA and the VCS decou- 
ple rather early and the derivations are rather lengthy. 
Below we provide a derivation of the relevant formulae 
starting from the output of a spectrometer. This allows 
us to present both the VGA and the VGS within a unified 
approach. Moreover, it allows us to cover the observa- 
tionally important case of crossing lines of sight. How- 
ever, unlike LP06, we do not consider self-absorption, 
which means that our approach is applicable either to 
optically thin emission lines or to weak absorption lines. 



3.1. Emission line 
Let us write out space-velocity distribution of density: 



no{r,v) = no{r)(p{v - Vmacroir)) 
where ^p(v) is Maxwellian distribution 

1 ^ ^ keT 



(1) 



ip{v) 



: exp 



(2) 



and nQ(r) is an integral density. 
If we write out the similar expression for emissivity 



e^{r,v) = eo(r)^e(u - 



ro(f)) 



(3) 



we will get Si^if) = ^ino{r), ipe{v) — ^fiiv) - for linear 
emissivity law, eo{r) = 72r7,Q(r), v'e(^) = 2^/t^lp'^{v) - 
for quadratic emissivity law, where 71, 72 set up the ratio 
between the quantities integrated over the velocity. 

Then we can write out the following expression for the 
signal in a spectrometer channel centered at the velocity 

f Wbi< 



SM = J ^^^^dr-£o(r),/(«.(r)+«r«(f)-z;o) 

(4) 



where 



/(«)= / ^,{v' + v)Uv')dv' 



(5) 



In Eq. [5] fs{v) is the channel sensitivity function, 
wi,{4i,9) is the instrument beam, eo(^) is the integral 
emissivity, Vr^r) is the turbulent velocity radius- vector 
projection, v^'^^{r) is the regular velocity radius- vector 
projection and A is the wavelength of the spectral line. 

In Eq. [4] it is assumed that the following normaliza- 
tion is applicable: J Wb{(f>, 9) dil = 1 and J_ fs{v)dv — 
1. °° 

The principal quantity in our analysis, S is measured 
in K, while Eq is in W/m^, A is in m and v is in m/s. 
We also assume that the instrument aperture efficiency 
is 1. 

Random field So{r) is not homogeneous, at least be- 
cause the emitting structure is limited in space. To model 
this we introduce a homogeneous field e{f) and a deter- 
ministic factor We{r) setting up borders of an observed 
object. We will also "pack" into e(r) all constant factors: 



£o{rj = ^3 We(r)e{r). 



(6) 



It is also convenient to introduce a window function w 
as follows: 



1 



w{r) = -^Wb{4',(>)we{r). 
With these notations we finally obtain: 



(7) 



= J w{r} dre{r)f{vr{r^ + vl''^{r) ~ v^) (8) 
In what follows we shall consider the Fourier transform 
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of spectral line^: 



S{K) 



1 

2^ 



S{vo)e' 



dvo 



= f{kv) J w{r)dfe{r)exp{-ik^{vr{r)+v^''^{r))) 

(9) 

This function can be easily determined from observa- 
tional data. 

3.2. Absorption line 

A signal in a channel dv received from a point object 
with the flux F is as follows: 



dP 



2kB 

2kB 



exp 



no (2:, v) dz ) dv 



1 - a y no(f)(y9(w - Vmacro{r)) dz j dv 

(10) 

Then the signal can be written as an integral of dP 
weighted with the channel sensitivity function fa : 

AeffF 



2kB 

l~a I no{z)f{v{z) + (z) - 1^0) dz 



Sa{vQ) 



And, if we make the following substitution^ 
AeffFa 



(11) 



2kB 



-no{z) = Wn{z)n{z) 



and keep only the variable term of So with inverted sign, 
we will have: 

5(«o) = j Wr,{z)dzn{z)J{vr{7^ + vl'S{r)-vo) (13) 
The correspondent Fourier transform is: 



S{ky) = f{ky) J Wn{z)dzn{z)ex.p{-iky{v{z)+v''''^z))) 

(14) 

The observational configuration for this regime is shown 
on Fig. H 

4. COMMON STATISTICS FOR VCA AND VCS 
Let us consider the following measure: 

Ki2{K,) = (Si{k^)S2 (ky) 

= f'^{ky) J wi{r)dr j W2{r') dr' 

[e{r)£{r')^ (exp{-iky{vr{f) - Vr'{r'))) 
■exp{-ik4v;'^ir)-vl.f 

^ Variable ky plays here the role of in LPOO, being however 
different in dimension (fc^ = bky). We use it here to avoid compli- 
cations when 6 = 0. 

^ as before, n{z) is a homogeneous field and Wn{z) is a corre- 
spondent window function 



(15) 



where indexes 1 and 2 designate different beam direc- 
tions. We have assumed here that density and velocity 
are uncorrelated. The effects of such correlation have 
been studied analyti cally in LPOO and numeri cally in 
iLazarian et"aH ()200 l) and Esquivel et"an (|2003l ). They 
were shown to be marginal. 

The first averaging gives us an emissivity correlation 
function Ce(r — r'). Averaging of the exponent can be 
performed with assumption that the velocity field has 
Gaussian statistics: 

exp{-iky{vr{r) - Vr'{r'))) 



exp 



To proceed we assume that the beam separation as 
well as the beam width is small enough to neglect the 
difference between Vr and Vz (we consider z-axis to be a 
bisector of the angle between beams). We also assume 
that w^'^f (r) depends only on z and admits a linear ap- 
proximation: 

vnz) = biz~zo)+v:';^. (16) 

This leads us to 
Ki2{ky)^Piky) [wiir)df [ W2{r')dr'C,{r-r') 



■ exp 



^Dz{T-?)-ikyh{Z-z')], 



(17) 



where is a velocity structure tensor projection: 

Dz{r- ?) = Uvz{r) - Vz{r')f) , (18) 



assuming that the velocity field is homogeneous. 

Having substituted r — r' , we can write the following 
expression for K12: 



Ki2{K)^ f{ky) / U.i2(f)df 



a(r)exp(-^D,(r) 



ik,,bz 



where 



1/712(7^)= / 'Wi{r')w2{r' + f) dr' 



(19) 



(20) 



which will be further referred to as "geometric term" . 

The case wuir) = 5{R — Ro)w^^a{z) leads us to the 
VCA^, while making the two beams coincide, we turn to 
the VCS. 

5. POWER SPECTRUM ALONG THE VELOCITY 
COORDINATE 

Setting wi = W2 in Eq. [19| we can write: 

Pl{ky) = Kl2{ky)\yj^^yj^ 

= ?(fc.) [wa(r)dr 



(21) 



C£(f)exp 



1^2 

'-^Dz{r) - ikybz 



^ Parallel lines of sight geometry is assum ed h ere, ij is a 2-d 
vector in picture plane, w^^a is defined in Sect l6.ll 
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Wa{r) = Wi2{r)\wi-- 



w{r')w{r' +r)dr' (22) 



Some variants of Wa are presented in Appendix [SI 
We get an important case of infinite resolution after 
applying of 5- function in Wa from Appendix 16.11 



(23) 

Ce(z) exp ( — ^Dz{z) — ikybz 

In particular, it arises when using an absorption line 
|13] from distant point sources behind the studied tur- 
bulent cloud (see Fig. 

Another important case arises when the emissivity 
spectrum is shallow and Cc splits into constant and sin- 
gular terms: 

(24) 



C,(f)«(£)^(l+(- 



and Pi is splitted accordingly: 

A(fc.) = Pl,„(fc.)+/'l,e(fc.) (25) 

In Eq. [25] the first term with the exception of its 
amplitude is entirely determined by velocity field, while 
the other one depends on both emissivity and velocity. 

6. GEOMETRIC TERM 

In this section wc consider geometric terms needed 
for specific variants of Pi. The cases presented be- 
low roughly cover all possible geometrical configurations, 
which include on the instrument beam shape and rela- 
tive location of the observed turbulent structure. Each 
of these variants of geometric term has different behav- 
ior for r ^ (see Tab. [5]), which results in different 
slopes for Pi asymptotics, even if velocity and emissivity 
statistics are the same (see Tab. [T]and Tab. [3]). 

6.1. Ideal resolution, pencil beam 

Let i? be a vector in a plane, perpendicular to z-axis. 
Assuming 6* « we have: 

Wb ^0, — ^ dR — 

So we can model "infinitely narrow" beam as 

Whir) = z^5{R) (26) 
Then the window function is as follows: 

w{r) K, \wb{r)w^{z) = S{R)we{z) 

and we can write the expression for Wa- 

Wa{r) = 5{R)We.a{z) 

where 

oo 

We,a{z) ^ / We{z')We{z' + z)dz' 



(27) 



(28) 



(29) 



Besides "very thin" pencil beam, this case is relevant 
for observation of turbulence in an absorption line of a 
remote point source (see Fig. [2]). 



6.2. Limited resolution, parallel lines of sight 

In this case emitting structure is far enough, so that the 
lines of sight do not converge within beam size. Therefore 
the corresponding window function is non-zero only 
in some vicinity of a distant point zq. We also assume, 
that picture-plane extent of the object is bigger enough 

than the beam size, so the dependence on R of can 
be ignored. 
Then we have: 

w{f) w \wb{9)we{z) 
4 

If we introduce 

Wb{R) = (- 
window function becomes 



/Wfc(i?)dJS=l (30) 



w{r) = Wb{K)we{z) 
and for Wa we can write: 

where 

Wb,a{R) = I Wb{R')Wb{R' + R) dR' 



(31) 



(32) 



6.3. Limited resolution, crossing lines of sight 

Here we consider the observation point to be inside or 
near the emitting structure. We also assume, that 
depends only on z. 

For our estimations we have to choose some particular 
form of Wb- For calculation convenience we may take the 
Gaussian: 



If we introduce 

VFb(f) E 




(33) 



z > (34) 
z < 



the window function takes the following form: 

w{r) = Wb{r)w^{z) 
and after some algebra we have: 

oc 

dz' 



Wa{r) 



1 



■ exp 



Z'2 + (Z'+|Z|)2 



i?2 



(35) 



92(2'2 + (z'+|z|)2)^ 
■We{z')w,{z' +\z\) 

6.3.1. Asymptotics 

As this function is singular at r = 0, we also need its 
asymptotic behavior. If Zedge is a characteristic scale of 
Wp, we can assume 



^OZedge 



and set — 
Wa(0, z) 



1. Then, if we set i? = 0, we have: 
1 



AOlz 



(36) 
(37) 
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Alternatively, setting z = we have: 

Wa{R,0) ^ 



2V2tt0oR 



(38) 



Approximating ^-dependence with an ellipse, we finally 
obtain: 



Wa{r) 



1 



1 



(39) 



2eor J 2718111^9 + 4:61 cos^ 61 



6.3.2. Approximation corresponding to tophat window 
function 

Equation [35| can be analytically estimated, if we take 
a tophat approximation for w^: 



w, 



_ / i> ze [zo,Zl] 
- 1 0, zi[zo,zi] 



(40) 



where zq and zi define the borders of the observed struc- 
ture. It gives us the following expression for the geomet- 
ric term: 



Wa{r) 



— 1 arctan ( 1 -|- 



2X0 ^ 



arctan ( 1 



221 



2^0oRz (2zl+pz^)~^ - {2zl-pz'^y^ 
R \ J R 



■ erf 



eo^2zi+pz-^ 



erf 



eo^2zi-pz'^ 



where 



P 



Zl + Zq 
Zl - Zo 



This case is applicable if the observation point is out- 
side, but still near to the emitting structure. This way 
we can account for the Local Bubble when observing the 
nearby turbulence in the Milky Way. 

6.4. Flat beam, parallel lines of sight 
Flat beam can be expressed as follows: 

Wb{e,x) = zwb,e{0)5{x) (42) 

In the case of remote emitting structure the corre- 
sponding window function is non-zero only in some 
vicinity of a distant enough point zq. Then the window 
function w is as follows: 

w{r) « —Wb,e{ — )5{x)w,{z) = Wb{R)we{z) (43) 

Zo Zq 

Then 



WbAR) = ^ / Wb,e{d)mA^+-)dO = S{x)Wb,eAy) 
Zo J Zo 



and 



Wair) = Six)Wb,g,aiy)We,aiz) 



(44) 
(45) 



6.5. Flat beam, crossing lines of sight 
In this case the window function has the following 



form: 



w{r) w -Wb,e{-)S{x)we{z) 
z z 



If we set 



wb,e{0) 



e 



and 



Weiz) = < e "■"'so , z > 
0, z < 



(46) 



(47) 



(48) 



we will have 



Wa[r 



6ix) 



dz' 



V^OoJ Jz'^ + (z' + |z|)2 





• exp 



+ {z' + \z\y 



9l{z'^ + iz' + \z\)^) 



edge 



(49) 



If we set y = 0, for z w we have: 

/ \ ^i^) 1 ^ Z 

Wa[X,0,z) « ^=;^l0g 



VZtT^O \Zedge 

Otherwise, if z = 0, for y « we get 
Wa[X, 0, z) w =— log 



Approximating dependence on 9 with an ellipse, we 
get: 



Six) 



2tt9o 



log 



^edge 



1 • 2 



COS^ 



(50) 



7. ASYMPTOTICS OF Pi 



The parameters of high-fct, asymptotics of Pi for dif- 
ferent emissivity spectrum types and geometries are pre- 
sented in Tab. ffl Tab. H TabU] and TabU Here "high" 
actually means "big enough for taking asymptotics for 
all remaining terms in Eq. [H]". Velocity structure ten- 
sor projection Dz in the exponent term is replaced by 
its asymptotics Eq. |A6| as well. Because of suppres- 
sion by this term, other terms are either factored out, if 
they are continuous at zero, or remain in the integrand, 
if singular. 

In particular, emissivity correlation function is con- 
tinuous for the steep spectrum and singular for the shal- 
low one (see Appendix IB.ip . Geometric term Wa is singu- 
lar in the regimes of infinite resolution (2-d (5-function), 
crossing lines of sight {r~^), flat beam with parallel lines 
of sight (1-d (5- function) , flat beam with crossing lines of 
sight (1-d (5-function and logarithmic singularity), and 
continuous for parallel lines of sight in low-resolution 
mode (see Sect. In general we can say, that for more 
prominent singularity in the integrand we have shallower 
Pi. 

The asymptotic regime of Pi depends on ky : with the 
real finite-size beam after some transition point we switch 
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from the high-resolution to the low-resolution mode with 
steeper slope. 

This point corresponds to the situation when the ex- 
ponent argument in Eq. |21| in transverse direction is 
about unity at beam radius zq0q. This gives us (for the 
parallel l.o.s.): 



For crossing l.o.s. with account for Eq. 
1 



It can be shown, that this is equivalent to 



(51) 

we have: 
(52) 



dn / k^dkul^{k) 



V 



(53) 



where rg = 2:0^0 for parallel and rg = Zedge(^ol'^ for cross- 
ing lines of sight. The latter expression is simply k^, 
which corresponds to the velocity variance at the beam 
scale. 

In Fig. [1] we illustrate the resolution regimes of the 
VCS technique. Depending on the beam size one may or 
not resolve the spatial extent of the eddies under study. 

To check the obtained results for asymptotical regimes 
we have performed direct calculation of the primary ex- 
pression of Pi [21]. The results presented on Fig. |4] 
confirm our estimations of asymptotic amplitudes and 
slopes. Numerical verifying with syntheti c field s is dis- 
cussed in detail in lChepurnov &: LazarianI ()2007D . 

8. DISCUSSION 

In the paper above we have presented the VCS tech- 
nique in a form that, first of all, reveals the close connec- 
tion between the VCS and the VCA. More importantly, 
we have described a new case very important for studies 
of galactic HI, namely, the case of collecting signal along 
converging lines of sight. 

We can change the effective instrument resolution to 
switch between regimes of Pi . This could be used as an 



additional test, if the chosen model is adequate. In prin- 
ciple, the relation between Pi amplitudes in high and 
low resolution regimes can give us additional informa- 
tion about ISM, such as compressibility. Further research 
should show how good must be the data to enable us to 
do such a study. 

It can happen that the asymptotic regime is not 
achieved because of thermal smoothing of the spectral 
line, or low spectrometer resolution. In this case we can 
calculate velocity field parameters by fitting of the model 
Pi to the data with different spacial smoothing. This 
appro ach is adopted, for instance, in IChepurnov et ahl 
(|200l . 

As mentioned above, VCS is a natural choice if we 
do not have an adequate spatial information. It can 
be applied to observations of absorption lines in spectra 
of a number of distant sources behind the cloud under 
study (Fig. [2]). Our simulations in a companion paper 
show, that the required number of independent measure- 
ments could be only ~ 10. The spectrum Pi is in high- 
resolution mode in this case, see Eq. [23] . In some cases 
we may be using the spectrum measured along only one 
line of sight, but use regular velocity shear to increase 
the effective statistical sample. In this case we "un- 
fold" the velocity fluctuations and this leads to reacher 
statistics, that we can utilize by averaging of adjacent 
fc^ -harmonics, see Fig. [3l 

Another important case for using the VCS can be the 
studie s of turbulence in galaxy clusters IVikhlinin et al.l 
(|1998f ). Such studies can utilize X-rays, but in most cases 
the spatial resolution of telescopes may not be adequate 
to get channel maps with a sufficient spatial resolution. 

As discussed earlier in LP06, the inertial range ob- 
served with the VCS may be somewhat limited. Some 
extension is possible by using heavier species, for which 
the thermal broadening is red uced. Another way that 
we explore in practical terms in IChepurnov et al.l (|2007l ) 
is fitting the data with a model and using the VCS for- 
mulae directly, rather than their power-law asymptotics. 
Naturally, the VCS is a complementary technique to the 
VCA and their use can improve the reliability of inter- 
preting fluctuations in the shapes of Dopplcr broadened 
lines in terms of underlying spectra of turbulence. 



APPENDIX 
VELOCITY STRUCTURE TENSOR 



In order to get the result normalized with respect to spectrum parameters we express the structure tensor projection 
Dz in terms of spectral tensor. 

If we co nsider v{r) so lenoidal with power-law power spectrum having cutoff at large scales, spectral tensor will be 
as follows (|Lesieuilll991| ): 



where Vq is a velocity power spectrum amplitude, a„ is a velocity spectral index and fcp is a cutoff wavevector. 
Then Dz is as follows: 

Dz{?) ^2jdk{l~- e'^^')zajF,j{k) 

After performing the integration over wavevector directions we will get: 
Dzir, 9) = 47rFoV"""3(/c(r) cos^ 9 + I,{r) sin^ 9) 



(Al) 



(A2) 



(A3) 



Velocity Coordinate Spectrum 
where Ic{r) and hir) are responsible for saturation of at large r's: 

dq ( (1 \2 1 2\ f , 3 / sing 

3 7 .^^"p^^^'^"^) ''^-v'vA~ 



W = i / 3rb exp Hkorf/q') . 1 - ^ ( ^ _ cos, ) ) (A4) 



OO 

/.W = 2/-i5L,xp(-MVfl-(i-^) 



CSO 

/-^expHvW).(l-^(^-cos, 



(A5) 



3 7 g"" 



If r 0, we can use 

D, (r, 61) w 47rVjfr"" (I^ cos^ 9 + 1, sin^ 9) ( A6) 

(in this text /» without argument means /*(0)) 
For potential field we have: 

F^^ik)-f^e (A7) 
In this case 



OO 

«'') = 2/^exp(-(...,W).(l-:i^) 



OC 



(AS) 



CSO 





EMISSIVITY CORRELATION FUNCTION 
Steep spectrum 

If spectral index > 3, the spectrum is long wave dominated and should have cutoff at low fc's to avoid divergence. 
If we assume an isotropic field, it may be as follows: 

F(fc) = ||cxp(-fcoVfc') (Bl) 
In this case the correlation function of a variable component of emissivity is 

OO 

C^Ar) ^ C,{^ - {sf = AnJk^dk^ exp {-kl/k') ^ (B2) 



The asymptotic behavior at r ~ is clear if we write as follows: 

Ce,o{r)^CM~ 47r£:oV"= -^Ist (r ) (B3) 

where 

OO 

dq ( ,2 2 I 2\ ( , sing 





For small r's we can replace Istir) with /st(0). 

Here we omit the second term in Eq. |B3| for the final asymptotics of Pi(fc„). 
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Shallow spectrum 

Short-wave dominated spectrum (a^ < 3) must have cutoff at large fc's: 

F(fc) = ||exp(-fcVfc?) (B5) 
The correspondent correlation function is 

oc 

C,,o(rO = AnJk'dk-^ exp (-fcVfc?) ^ (B6) 



It can be rewritten as follows: 
where 

oo 

/..W=/^expHVfc^^)^ (B7) 



Considering ki being high enough, and keeping in mind, that relevant r's should significantly exceed** , we can 
use the following approximation: 

C^Ar) « ^/.^ (B8) 

where 

oo 


Quadratic emissivity 

The expressions |B1| and |B5| give us realistic approximation for low fc's for linear emissivity. In the case of quadratic 
one self-convolution of spectra of such forms should be used instead: 

F,{k) = J Fp[k')Fp{k' - k) dk' (BIO) 

It is quite clear that for steep density the emissivity spectrum has asymptotically the same slope, than the underlying 
density one. However, for shallow density we have a different picture: 

F,{k) ~ J |^fc;2 + ^ j^fc, _ ^ ^ y2 ^ (^y^ ^ ^ 

oo 1 / 



2^ / dk' j [{ k'^ + -kA ' - k'^khA dt 









1 



(Bll) 



Therefore, if emissivity is proportional to squared density, emissivity spectral index is as follows: 



* These values of r can however be enough small to take r ^ asymptotics for other terms in Eq. I21| 
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ACCOUNTING FOR THE LINEAR VELOCITY SHEAR 

Let us consider the case of ideal resolution. With account for Eq. [55] and keeping the regular velocity term, we will 
have Eq. |21| in the following form: 

Pi{K)^?{K) J w,,a{z)dzC,{z)exp(^-'^D,{z)~iKbz^ (CI) 

— oo 

Or, for high ky^s^ 

P,,.j.H^i4<2)*..i!<M (C2) 

where G{ky) is a factor, responsible for the hnear velocity shear: 

oo 

G{K) « — / exp f-g""-^) cos dq (C3) 

where 

H = i2nV„^Ic)^b-~^ (C4) 
It is obvious, that 

lim G{ky) = 1 , 3 < a < 5 

so the high-fc^ asymptote stays valid in this case, too. 
The same can be easily shown for the case of limited resolution. 



CONSTANTS, INVOLVED IN Pi ASYMPTOTICS 



oo 

Jz,i = yexp(-g""-3) dq (Dl) 



oo 

^^^^-/tSt-pH--^) (D2) 







oo 

Jr,l= J exp{~q°'^'~^)q^dq (D3) 


TT 

f sin0d9 
{Ic cos'^ 9 + 




oo 



Jr,2 = yexp(-g"--3)g"--id(z (D5) 



} sm9d9 , , 

y (/cCos2 6l + /,sm26')^ 

Jr,3 = J exp{~q''^-^)qdq (D7) 


} sm9de , , 

J9,3 - / , — — (D8) 

{ J2tt sin^ 9 + 4612 cos2 0{I,sin^0 + 4 cos2 9)^^ 
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Jr,4 = J exp (-(j""'"^) 9"="^ 



(D9) 



Jr,[ 



{ J 271 sin^ e + 461^ cos2 9{Issin^9 + h cos^ 0)^- 



exp (— q"" g 



(DIG) 



(Dll) 



J9.I 



/ 



00 



d0 



{IgSin'^O + Ic cos^ ( 



J,,6 = y"exp(-g"^-3^ 



(D12) 



(D13) 



de 



{IsSiri^O + /c cos^ ( 
Jr,i = y ^'^P {-q""'^^) qdq 




00 



d6i 




00 



log 



Jr,8 = / exp {~q 



J- sin^ e + cos2 61 

Wo 



(DM) 
(D15) 

(D16) 
(D17) 



Je.i 



dO 



Ic cos^ I 



-log 



TT sin^ e + cos2 61 



edge{2Wi)—^{h cos^ 9 + h sin^ ^)^ |A:„| 



(D18) 



Here the first index designates the integration domain. 
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TABLE 1 

ASYMPTOTICS OF Pl(kv) AND Pl^vikv) FOR kv —* OO 



l.o.s. geometry pencil beam 



flat beam 



low resolution 



parallel 



crossmg 



-2 

lfc,.l 



4, ^ 



I 



ei 

4 



Note. 



Amplitudes Ay see in Tab. [2] 



TABLE 2 

Normalization coefficients for Pi(fc„) and Pi,i,(fc^) 



l.o.s. geometry pencil beam 



flat beam 



low resolution 



parallel 
crossing 



2me,a(0)<Tj: Jj_i 
2me,„(0)aeJz,l 



'y!,.£).a(0)^'.,a(0)'^eJr.5-^ 
■ — ■ 2 



(2^vo2)^j^;^=^ 

'2 



Note. — For Pi = Ce(0), for Pi,„ ctj = (e)^ 



TABLE 3 

Asymptotics of Pi,E(fc^) for fe„ ^ 00 



l.o.s. geometry 


pencil beam 




flat beam 


low resolution 


parallel 










Qe-2 




-1 




crossing 




|fc„| oiJ- 
A,p{k^ 




|fc„|°"-^ 






ifc„r°--3 





Note. — Amplitudes Ae see in Tab. |4] 



TABLE 4 

Normalization coefficients Ae for Pi,e(fc«) 



l.o.s. geometry pencil beam 



flat beam 



low resolution 



parallel 



8irtUe,a(0)£o-f5h-^2,2 2ir£o-fJh Jr.S Je.S 
7^7^ 



(2,rV(2)-^.-3 
(2^Vo2)-i,-3 



TABLE 5 

Asymptotics of geometric term Wa{r) for r ^ 



l.o.s. geometry pencil beam 



flat beam 



low resolution 



parallel 
crossing 



"'^,a(0)14^5,„(0) 

1 



Bar ^2n sin^ e+4ef^ cos^ B 



Note. — f= (x, y, z), R = {x, y), 9 is an angle between r and z. 
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TABLE 6 
List of symbols 



symbol 



description 



defined at 



k 
R 

V 

Ctv, O, 

no 

eo, e 
S 

fs,f 

X 

Wt 
We 
W 

K 

Ce 

W12, f 
Pi 

Ic, Is 
J.,. 



radius- vector 

Fourier counterpart of r 

picture-plain projection of r 

velocity 

Fourier counterpart of v 

spectral indices for velocity, emissivity and density 
number density 
Maxwellian distribution 

emissivity, real and idealized (homogeneous) 
spectral line 

channel sensitivity function, original and effective 
wavelength of a spectral line 
instrument beam 

factor, defining the observed object boundaries 
window function 

statistics, common for for VGA and VCS 
velocity structure tensor projection 
emissivity correlation function 

geometric term: general, for coinciding beams, for 
power spectrum along velocity coordinate 
auxiliary function, derived from Wi, 
auto-convolution of Wt 
auxiliary functions for Dz 
auxiliary constants 



Sect. |3] 
Appx. |A] 
Sect, m 
Sect. 
Sect. |3] 

AH . IBII 

1] 




Appx. EH 



M 

M< ED 

M 

A4] . [A5] 
Appx. IDI 



Note. — Tilde means Fourier tra nsfo rm, v and e subscripts in Pi and K mean velocity 
and emissivity terms according to Eq. 1251 . 



Eddie modes: 

1 - low resolution 

2 - high resolution 

3 - intermediate 



Cloud 



Beam 





V 



Fig. 1. — Illustration to VCS resolution regimes. Eddies within beam size are in low resolution mode and cause steepening of Pi. The 
ones exceeding the beam size are in high-resolution mode and correspondent Pi is shallower. Therefore high-resolution mode is more 
preferable regarding to the signal-to-noise ratio. 
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Fig. 2. — VCS allows recovery of velocity statistics from absorption lines from stars. The whole Pi is guaranteed to be in high-resolution 
mode, which provide better dynamical range over ky. Numerical simulations show, that very few independent measurements are needed 
to gain enough statistics. 
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Improving statistics by regular velocity shear 




-'/'^ -?.Ti -■.= -■ ^ -1.? 

One realization, spectaim compression is 16. 
Aif due to galaxy rotation is 8 times greater 
than turbulent line width. 



Fig. 3. — VCS can be applied even to a single spectral line. Configuration can be as shown on this picture. Regular velocity shear from 
the galactic rotation results in statistics being sufficient to recover Pi 




Fig. 4. — One-dimensional power spectrum along velocity coordinate (direct calculation). Left: parallel lines of sight, distance to the 
object is 1000, object thickness is 100, injection scale is 30. Right: observation point is inside the emitting structure (converging lines of 
sight), distance to emitting structure's edge is 1000, injection scale is 30. 



